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ABSTRACT 

We present the observations of the pulsar SMC X-1 with the ART-P telescope 
onboard the Granat observatory. We investigate the variability of the flux from the 
source on time scales of several tens of days. The intensity variation of the pulsar 
are shown to be consistent with the presence of a periodicity in the system with a 
characteristic time scale of ~ 61 days. The precession of an inclined accretion disk, as 
indirectly confirmed by the absence of low-state pulsations, may be responsible for the 
observed variability. The spectrum of the source is well described by a power-law energy 
dependence of the photon flux density with a slope of ~ 1.5 and an exponential cutoff at 
energies above ~ 14 — 18 keV. We estimated the inclinations between the planes of the 
orbit and the accretion disk and the magnetic field of the neutron star. 
©2004 MAIK "Nauka/Interperiodica". 
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INTRODUCTION 

The X-ray pulsar SMC X-1 located in one the nearest galaxies, the Small Magellanic 
Cloud (SMC), is one of the most intense and rapidly rotating accreting X-ray pulsars. 
The source, which was first discovered in X-rays in 1971 (Price et al. 1971), is a member 
of a high-mass binary together with the companion star Skl60, a BO-type supergiant with 
a mass of 17.2Mq and a radius of 18Rq (Reynolds et al. 1993). On short time scales, 
its flux pulsates with a period of about 0.71 s, while the neutron star itself is eclipsed by 
its companion for ~ 15 h (the duration of the X-ray eclipse) with a period Porb — 3.892 
days. Analyzing the results obtained over the entire preceding history of observations 
of this source, Levine et al. (1993) determined the orbital parameters of the binary and 
estabhshed that its orbital period is decreasing at a rate Porb/ Porb — —3.35 x 10~^ yr~^- 
Apart from the rotation period of the neutron star and the orbital period, the observations 
are indicative of a 50-60 day cycle in the binary (Wojdowski ct al. 1998), during which the 
X-ray luminosity of the source varies between 10^^ erg s"^ and several units of 10^^ erg 
s~^. The latter value is close to or even higher than the Eddington limit for a spherically 
accreting neutron star with a mass of IAMq. This fact as well as the observed steady 
spinup of the neutron star for the more than 30 years since the discovery of the pulsar (see, 
e.g., Wojdowski et al. 1998; Bildsten et al. 1997; Lutovinov et al. 1994; and references 
therein) suggest that a disk accretion model is most hkely reahzed in the Skl60/SMCX-1 
system; the accretion disk itself can be formed either from stellar wind or during mass 
transfer through the inner Lagrangian point. Analyzing the EXOSAT light curves for the 
pulsar SMC X-1, Angelini et al. (1991) found a sharp (by a factor of ~ 3), short-time 
increase in the flux from this source with a duration of ~ 80 s, which they interpreted as 
a type-11 X-ray burst. Recent RXTE observations of the pulsar have revealed more than 
a hundred such events with time scales of several tens of seconds (Moon et al. 2003), 
suggesting that the radiation and accretion mechanisms in the source SMC X-1 under 
study and the pulsar/burster GRO J174428 are similar. In this paper, we analyze the 
archival data obtained with the ART-P telescope onboard the Granat observatory to 
study the long-term flux variations of the source as well as its temporal and spectral 
parameters. 

OBSERVATIONS 

The ART-P X-ray telescope onboard the Granat international orbital astrophysical 
observatory consists of four coaxial, completely independent units, each of which includes 
a position-sensitive detector with a geometric area of 625 cm^ and an URA-based coded 
mask. The telescope can image a selected region of the sky by a coded-aperture technique 
in a 3°. 4 x 3°. 6 field of view (the full beam width) with a nominal resolution of ~ 5 arcmin 
in the energy range 3-60 keV. It has an energy resolution of ~ 22% in the calibration 
5.9-keV iron line. The observations were carried out in the photon-flux mode. In this 
mode, for each photon, its coordinates on the detector, the energy (1024 channels), and 
the arrival time (the accuracy of the photon arrival time is 3.9 ms, and the dead time is 
580 s) were written to buffer memory. This mode allows one to perform both timing and 



Table 1. ART-P observations of the pulsar SMC X-1 in 1990-1992 . 
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After correction for the motion of the Solar-System barycenter and the orbital motion in the binary. 
In the energy range 620 keV. 

In the energy range 620 keV for the assumed distance of d = 50 kpc to the source. 
For technical reasons, the spectrum of the source cannot be restored. 
The 3cr upper limit for the pulse fraction is given in mCrab. 
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spectral analyses of the radiation from each X- ray source within the field of view of the 
telescope. A detailed description of the telescope was given by Sunyaev et al. (1990). 

The Granat observatory observed the region of the SMC that contained the pulsar 
SMC X-1 in series once a year (Table 1). A total of three series of observations with a 
total exposure time of ~ 271000 s were carried out for the source, which allowed us to 
study in detail the radiation from the pulsar as well as its spectrum and variability on 
various time scales. Preliminary results of our timing analysis of the ART- P/Granat 
observations of the pulsar SMC X-1 were presented previously (Lutovinov et al. 1994). 
The typical pulse profile for the source in the energy ranges 3-6, 6-10, 10-20, and 20-30 
keV (Fig. 1) exhibits two symmetric peaks located in the phase ranges 0.1-0.4 and 0.5-1.0, 
with the width of the second peak decreasing with increasing photon energy. 

It should be noted that during the first series of observations (in the spring of 1990), 
we used the first module of the telescope. The subsequent observations were performed 
with the third module, which had a lower sensitivity at soft energies; therefore, all of the 
results presented here refer to the energy range 6-20 keV. In addition, the high-voltage 
level of the detector of this module varied significantly over time, and it was not always 
possible to carry out calibration observations of the Grab Nebula. As a result, we ran 
into considerable difficulties in constructing the response matrix for some of the observing 
sessions, making it difficult or, in several cases, even impossible to carry out a spectral 
analysis of the X-ray radiation recorded with this module. 

LONG-TERM INTENSITY VARIATIONS OF THE SOURCE AND 

BURSTS 

As was noted above, the long-term observations of SMC X-1 showed that, apart from 
the proper rotation of the neutron star and its orbital motion, the system has another, 
nearly periodic component-intensity variability of the source on a time scale of several 
days, with the recorded X-ray flux decreasing by more than an order of magnitude. 
A similar pattern is also observed in the other two well-known X-ray pulsars Her X-1 
and LMC X-4, where it is attributable to periodic eclipses of the emitting regions of the 
neutron star by a precessing accretion disk (Lutovinov et al. 2000; La Barbera et al. 2001; 
and references therein) . It is assumed that the same mechanism may also be responsible 
for the observed intensity variations in the pulsar SMC X-1. 

To determine the precession period Pprec of the accretion disk, we fitted the AKT-P 
data by a sinusoidal signal with a trial period varying in the range of 40-80 days. The 
deviations of the measured fluxes from their predicted values were determined by the least- 
squares method. The best value was obtained for a period Pprec — 61 days. This period 
agrees with the results by Wojdowski et al. (1998), who found long-term variations in 
the X-ray flux from SMC X-1 with a period of 50-60 days when simultaneoiisly analyzing 
the ASM/RXTE and BATSE/CGRO data. The limited set of data (observing sessions) 
makes it impossible to completely cover the entire presumed period. Therefore, our 
period estimate is not an accurate and statistically significant measurement of the possible 
precession period Pprec (the statistical significance of the peak on the periodogram is 
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~ 1.5a). It most likely gives circumstantial evidence for the presence of a third type of 
periodicity in the system and the possible mechanism of its formation (see below). 

The analysis of the light curve for the source by Clarkson et al. (2003) shows that the 
presumed precession period is not constant in itself, but varies smoothly over an interval 
of 40-60 days with a characteristic time scale of ~ 1600 days. 

Figure 2 shows the light curve for SMC X-1 constructed over two years of ART-P 
observations in the energy range 6-20 keV. The dots with (la) error bars indicate the 
measured fluxes from the source in mCrab during individual observing sessions, and the 
solid line represents their best sinusoidal fit with a period of ~ 61 days. It is undoubtedly 
of considerable interest to compare the phases of the light curves obtained simultaneously 
with the ART-P telescope and the BATSE observatory. However, the quality of the latter 
in this period was too low to make such a comparison (sec Fig.l from Wojdowski et al. 
1998). A comparison with more recent RXTE observations of the pulsar has revealed 
no correlation with our results, which is most likely attributable to the variability of the 
precession period (see above). 

When analyzing the light curves for the presence of X-ray bursts, we found several 
events that were similar to those recorded by the RXTE observatory (Moon et al. 2003) 
and that could be interpreted as weak type-II bursts. Note that the statistical significance 
of the recorded events is low, 3—5(7. However, their roughly triangular shapes and HWHM 
of ~ 13 s obtained by Gaussian fitting are similar to those observed previously. Figure 
3 shows the profile of one of such bursts. 



The spectrum of the X-ray pulsar SMC X-1 is typical of this class of objects and can be 
described by a simple power law. In general, depending on the specific source, this model 
is modified by an exponential cutoff at high energies, absorption at soft energies, and 
emission or absorption lines. In the case of SMC X-1, we did not detect any significant 
features in emission or absorption. In addition, as our study shows, the interstellar 
absorption estimated by different authors to be in a wide range, 5.9 x 10^° — 3.2 x 
10^^ atoms cm~^ (Wojdowski et al. 1998; Moon et al. 2003), plays no important role in 
our analysis of the radiation from the X-ray source in the energy range 6-20 keV. Thus, 
the ultimate formula that we used in our spectral analysis of the radiation from SMC X-1 
is 



Here, E is the photon energy in keV, Iiq is the normalization of the power-law 
component to 10 keV, is the photon spectral index, Ec is the cutoff energy, and Ef is 
the characteristic e-folding energy in the spectrum of the source (although this is a purely 
empirical formula with a cutoff at Ec, it is widely used to fit the spectra of X-ray pulsars; 
see White et al. (1983)). 

Figure 4 shows the photon spectra for the pulsar (in photons cm~^ s~^ keV"^) 
measured during several observing sessions with different intensity levels. The curves 
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represent the best-fit model spectra for tfie source based on a simple power law or formula 
(1); the best-fit parameters are given Table 2. In choosing between a simple power law 
and its modification by a cutoff at high energies, we used the A^^-test. Based on this 
test, we determined the probability that it was no accident that the value improved 
when passing to a more complex model. The statistical significance of this passage is 
higher than 95% in most cases. 

If we assume that the low state is observed when the source is hidden behind the 
outer edge of a warped or inclined (to the orbital plane) accretion disk and that we 
see the flux attenuated by absorption and scattering in the cold material on the disk 
periphery or scattered in the hot corona above the disk (as is the case in the system Her 
X-1 (Lutovinov et al. 2000)), then, to a first approximation, the spectra of the high and 
low states may be considered to differ only by normalization and additional absorption. 
Fitting the low-state spectrum of the pulsar (April 5, 1990) by formula (1) with the 
parameters fixed at their high- state values does not give a positive answer to the above 
assumption (an appreciable increase in ). However, it cannot be completely rejected 
either, because the data are statistically limited. In this case, an important criterion 
could be an increase in absorption at low energies. However, we cannot be certain about 
this, because the energy range is limited: the measured absorption in the low state is 
estimated to be (31.5 ± 12.0) x 10^^ atoms cm~^, while the upper limit for the absorption 
is 20 X 10^^ atoms cm"^ (1<^)- 

To investigate the spectral evolution of the source on a scale of one pulse period, we 
performed phase- resolved spectroscopy of the pulsar radiation in the high state for two 
observing sessions, April 24 and May 11, 1990. To this end, the data obtained during 
each observing session were folded with the pulsation period and divided into 12 time 
bins (the arrival time of each photon was first corrected for the barycenter of the Solar 
System and the motion of the neutron star in the binary). The total accumulation time 
of the signal for each of the 12 spectra was 1137 and 1568 s for the sessions of of April 24 
and May 11, respectively. The subsequent analysis indicated that the spectra obtained 
in this way are well described by a simple power law in the energy range 6-20 keV; the 
spectral slope is virtually independent of the pulse phase. This conclusion is illustrated 
by Fig. 5, which shows a plot of the photon spectral index against the pulse phase for 
the observing session of April 24. 

DISCUSSION 

High and Low States 

It follows from the observations of SMC X-1 by different observatories that the 
recorded X-ray fiux from this pulsar is not constant but undergoes large, nearly peri- 
odic variations on a time scale of 50- 60 days. By analogy with the well-known binaries 
SS 433, Her X-1, and LMC X-4, which exhibit a similar pattern, the precession of an 
inclined accretion disk is considered to be one of the main causes of the observed vari- 
ability. The normal companion Skl60 in the binary under consideration is a supergiant 
with an intense stellar wind that provides the bulk of the accreting material forming 
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the accretion disk. The mechanism of its formation is not yet completely understood, 
because, in general, supergiants are close to filling their Roche lobes and material can be 
transferred to the compact object through the inner Lagrangian point. Both mechanisms 
may operate in the system simultaneously. It should be noted, however, that, according 
to Clarkson et al. (2003), even if the stellar wind plays a major role in the formation of 
the disk, it is a coUimated wind, and, thus, the two cases are characterized by similar 
patterns of mass transfer to the accretion disk. 

Larwood (1998) showed that a relation exists between the precession period and the 
binary parameters for a precessing accretion disk: 

^ = (3/7)g(l + q)-'/\Ro/af'\os5, (2) 

-'prec 

where q is the mass ratio of the normal component and the compact object; Ro is 
the outer radius of the accretion disk, which may be expressed in fractions (5 of the 
Roche lobe size; a is the separation between the binary components; and 5 is the angle 
between the orbital and disk planes. It follows from an analysis of the optical light 
curves for Skl60/SMC X-1 (Howarth 1982; Khruzina and Cherepashchuk 1987) that the 
size of the accretion disk is /9 ~ 0.7 — 1.0. For the mass ratio q = 10.8 and periods 
Porb/Pprec — 0.064, wc cau estimate the range of possible inclinations, () 25" — 58°. 
The minimum angle between the orbital and disk planes, cos6 ~ 1, corresponds to an 
accretion-disk size of /3 ~ 0.65, which slightly exceeds the maximum disk radius, P ~ 0.61, 
obtained by Paczynski (1977) for accretion only through the inner Lagrangian point. This 
discrepancy is even larger if we take into account the range of possible precession periods. 
Thus, mass transfer from the normal companion only through the inner Lagrangian point 
in the system Skl60/SMC X-1 is unlikely. 

As was noted above, no X-ray pulsations were found during the low state of the source, 
and only upper limits for the pulse fraction were obtained (Table 1). However, since these 
limits are rather large, we cannot completely rule out the presence of pulsations in the 
low state. Wojdowski et al. (1998) presented ROSAT measurements of the pulsation 
period in the low state. However, they stipulate that these measurements were made at 
the very beginning of the low state, when the source may have not yet been completely 
shielded by the disk. A similar situation was probably observed by the ART-P telescope 
on May 18, 1992, when the intensity of the source was low, but X-ray pulsations were 
recorded with low statistical significance (~ 4cr). 

Magnetic Field 

Based on the observed parameters of SMC X-1 obtained over a long observing period, 
we can try to estimate the magnetic field of the neutron star in the system by using the 
model of an accretion disk suggested by Li and Wang (1996). 

The history of measurements of the pulsation period for SMC X-1 is indicative of a 
virtually uniform spinup of the neutron star since its discovery. The mean rate of change 
of the period during our observations in the period 1990-1992 was P = —(32.6 ± 0.8) x 
10"^ s yr~^. Subsequent ROSAT, ASCA, and RXTE observations of the source yield 
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a similar value, P = -(32.0 ± 0.2) x 10"^ s yr"^ (Wojdowski et al. 1998). Thus, 
assuming that the observed variability of the source is not related to the variations in the 
intrinsic radiation from the pulsar, but is determined by extraneous effects, we estimated 
its bolometric luminosity from the observed spectral parameters as L^; = 47 x lO^'' erg 

In a situation where the neutron star spins up almost uniformly over a long time 
interval, the change in its rotation period is related to the binary parameters by the 
equation 

._ M(GMro)V^n(^,)p^ 

where M is the accretion rate; / is the moment of inertia of the star; tq is the 
inner edge of the accretion disk, which in our model is assumed to be equal to the 
Alfven radius r^; n{ujs) — 1 + ^^31^1^^"^)"^ is the dimensionless angular momentum; and 

— ^s/^fc('"o) is the speed parameter. The latter depends on the angular velocities 
Qs and Jlfc(ro) of the neutron star and the material at the inner edge of the accretion disk, 
respectively. We estimated the moment of inertia of the neutron star from the relation 
/ = OAM^R"^ and assumed that its radius is i? = 10^ cm, the mass is = 1.4Mq. 

The magnetic moment fi of the neutron star estimated using formula (3) shows that 
the best agreement with the experimental data is achieved for n ~ (0.05 — 0.1) x 10^° G 
cm^, which is equivalent to a neutron-star surface magnetic field of ~ (1 — 2) x 10^^ G. 
This conclusion about the relative weakness of the magnetic field is consistent with the 
observation of bursts from the source that can be classified as type-II bursts (see Moon 
et al. 2003). Moon et al. (2003) also estabhshed that these bursts have much in common 
with the bursts observed previously from another super-Eddington pulsar, GRO J1744- 
28. They suggested separating out the so-called group of pulsars/bursters with magnetic 
fields of 10^^ G. Another important argument for this hypothesis is the absence of 
cyclotron features in the spectrum of the source in the energy range 3-100 keV, which 
corresponds to magnetic fields of > 5 x 10^^ G. 

Note that the applicability of the model under consideration at large is limited by 
the fact that the speed parameter reaches the critical value cUc at which the angular 
momentum is not transferred to the neutron star. Having reached this level, the pulsar 
must pass from spinup to spindown, which is not observed for SMC X-1. The fact that 
approaches 10^° G cm^ means that the speed parameter falls within the range of critical 
values cus- 



CONCLUSIONS 

The X-ray pulsar SMC X-1 was repeatedly observed in 1990-1992 with the ART-P 
telescope onboard the Granat observatory. Over this period, we accumulated the data 
that allowed us to investigate the variability of the source on various time scales and its 
spectrum and to estimate the binary parameters. 

We showed that, apart from the periodicities associated with the proper rotation of 
the neutron star and the orbital motion, the system has yet another, nearly periodic 
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componcnt that may be associated with the precession of the accretion disk. Its period 
Pprec ~ 61 days, as inferred from the ART-P data, agrees well with the observations of 
other observatories. 

The spectrum of the source is typical of X-ray pulsars and can be described by a 
simple power law with a cutoff at high energies. Our analysis showed that the spectral 
shape depends weakly on the orbital phase and intensity of the source. Phase-resolved 
spectroscopy of the radiation from the source in the high state revealed no dependence 
of the spectral slope on the pulse phase either. 

The above estimates of the magnetic field for the neutron star show that its strength 
must be ~ 10^^ G to be consistent with the observational data (the spinup rate the star. 
X-ray bursts, and the absence of cyclotron features in the spectrum). 
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Fig. 1: Energy dependence of the pulse profile for SMC X-1, as derived from the ART-P data 
of April 24, 1990. The errors correspond to one standard deviation. 
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Fig. 2: The 6-20-keV light curve for SMC X-1 constructed over two years of ART-P observa- 
tions. The dots indicate the measured fluxes during an individual observation, and the solid 
line represents the best sinusoidal fit with a 61-day period. 
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Fig. 3: The profile of a typical X-ray burst detected by ART-P from SMC X-1 during the 
observing session of April 22, 1990. Time in seconds from the beginning of the observing 
session is along the horizontal axis. The time bin is 8 s. 
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Fig. 4: The photon spectra for SMC X-1 measured with the ART-P telescope during the first 
series of observations (the spring of 1990). Different symbols (circles, squares, and triangles) 
indicate the photon spectra for the source in states with different intensity levels; the curves 
representtheir best fits by a power- law decrease in photon flux density or formula (1). 
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Fig. 5: Photon spectral index versus pulse phase. The errors correspond to one standard 
deviation (a). 



